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The activation energy of viscoelastic deformation of rocks 


DZHEMA V. MELKONYAN 


Abstract. The suggested model, which considers the simultaneous appearance (on the intermolecular scale) 
of elastic and viscous behaviours of rocks during their viscoelastic deformation, enables the relaxation time, tọ 
of the viscoelastic deformation energy to be defined. By determination of tọ, the activation energies of the elas- 
tic and viscous deformation can be separated from the temperature dependence of the viscoelastic deformation 
energy stream U. The probabilistic approach to viscoelastic deformation of rocks allows the introduction of the 
a notion of the total probability, consisting of the probabilities of elastic and viscous deformations. 
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Абстракт. Модел који je предложен у овом раду приказује симултану појаву (у интермолекуларној pa- 
змери) еластичног и вискозног понашања стена током њиховог вискоеластичног деформисања. На тај на- 
чин омогућено је дефинисање времена релаксације tg вискоеластичне деформационе енергије. Када се ово 
време одреди онда могу да се раздвоје активиране енергије еластичних и вискозних деформација од тем- 
пературне зависности вискоеластичног деформационог енергетског тока И. Пробабилистички приступ, 
примењен на вискоеластичне деформације стена, омогућава да се уведе појам тоталне вероватноће која 
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се састоји од вероватноћа еластичних и вискозних деформација. 


Кључне речи: вискоеластичност, деформација, активирана енергија, енергетски ток, време релаксације. 


Introduction 


At present, there is great interest in the viscoelastic 
deformation of rocks, since the viscoelastic bodies are 
heterogeneous and complex, in particular the deforma- 
tion of rocks which elapse at extraordinarily various 
physical conditions existing in the interior of the earth 
(e.g. TROITZKAYA, 1961). 

The crust rocks at different natural conditions are 
both elastic and brittle solids. The ability of crust and 
the mantle rocks to transmit seismic waves, when the 
forces act during very short times, testifies to their high 
elasticity of compression and displacement (e.g. VER- 
HOOGEN еї al., 1974, Товсотт & SCHUBERT, 1985). 
However, there is evidence for the elastic behaviour of 
rocks at more long stresses, for example Chandler oscil- 
lations of the pole and earth tide. On the other hand, 
the brittle behaviour of rocks in the crust upper is also 
verified by the existence of numerous natural break sur- 
faces on all scales (breaks, cracks, etc.), (e.g. WERHO- 
OGEN et al., 1974, TURCOTT & SCHUBERT, 1985). There 
are other data (e.g. VERHOOGEN et al., 1974), which 


confirm: the viscous flow of rocks in the solid state 
without large breaks; the frequent appearance of inelas- 
tic and residual natural deformations in rocks. The ten- 
dency of the crust to a hydrostatic equilibrium allows 
the assumption that some uppers of the earth’s crust be- 
have like a very viscous liquid in their properties. Crust 
rock flow in the solid state on a small scale is proven 
by such structures as foliation, orogeny and the defor- 
mation of bodies with certain primary sizes or shapes. 

The above stated data point to the fact that elastic 
and viscous deformations occur simultaneously in the 
crust. 

In the theory of macroscopic elasticity of viscoelas- 
tic deformation, the value of the characteristic time is 
usually given, during which the elastic behaviour changes 
into viscous behaviour. This is a rough model of rocks: 
elastic under the influence of momentary forces and 
viscous under the influence of the long duration of tec- 
tonic forces. In reality both behaviours (elastic and vis- 
cous) appear simultaneously (e.g. ZELINSKY & MELKO- 
NYAN, 1997). TROITZKAYA (1961) and VERHOOGEN et al. 
(1974), assumed that the elastic behaviour of a body is 
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slowed down by the viscous component. For example, as- 
sume that during a given time (let it be assumed, ~ 100 
years), the elastic behaviour of a substance completely 
changes into viscous. However the viscous behaviour of 
a deformable body appears during the deformation pro- 
cess (on the intermolecular scale), (e.g. LANDAU & LIF- 
SHITS, 1987), though the body is elastic as a whole. Such 
an assumption allows the fact that the intermolecular 
forces have an insignificant radius of action (e.g. LANDAU 
& Lirsuits, 1987). In spite of the fact that the influence 
of the forces extends around, thus making them particles 
only at an intermolecular distance, it simultaneously cre- 
ates a residual phenomena, which can be involved even 
at the beginning of the deformation processes. 


The energy stream accumulation and 
relaxation at viscoelastic deformation 


The purpose of this work was to investigate the vis- 
coelastic deformations of rocks and to determine their 
activation energies taking into account the simultanious 
apperance of elastic and viscous behaviours and defor- 
mation energies of rocks. The viscoelastic deformation 
is considered in certain deformable plane of rocks with 
thickness /. 

For a deformation €, (e.g. VERHOOGEN et al., 1974), 
the following equation can be written: 


e(x,t) _ 1 Or(x,t) e(x,t) 
а п Ox to 


а) 


where n is the viscosity, т is the stress, tọ = n/u is the 
relaxation time of viscoelastic deformation, where is 
the shear modulus. From the time tọ, the relaxation time 
of the energy of viscoelastic deformation can be identi- 
fied. The axis x is normal to the deformable plane. 

Eqn. (1) is valid for linear deformation, i.e. tọ = const. 
Taking into account the deformation energy, Eqn. (1) 
can be integrated with respect to the thickness / of the 
deformable domain from which the cumulative elastic 
energy is intend to be released, and by multiplying to 
the area S of this domain, one obtains an equation de- 
scribing the energy stream accumulation and the relax- 
ation in time in the deformable domain: 

ао о 


a> eee (2) 


in which the energy stream on the boundary of the de- 
formable domain and the bulk particles is given by: 


A 
ООВ (3) 
Ш] 


where 7(0,f) denotes the total stress on the plane x = 0, 
ће. on the mentioned boundary. Taking into account 
that the influence of the shear resistance on the accu- 


mulation and relaxation of the elastic energy is negli- 
gible, ie. т, = 0, then A,(t) will iterate the shape of 
the elastic stress through the deformable domain. For 
an energy stream which is released from a deformable 
domain, one can write: 


д 
и = on Q (4) 
ЈЕ 10 
in which y; denotes the carry coefficient of the ener- 
gy, which can be expressed as: 
ТЕ 


УЕ = (5) 


ТЕ +Ти 


where Тр and т, are the elastic and viscous stresses, 
respectively. The physical meaning of the coefficient 6, 
will be explained below. 
Substitution of expression (4) into Eqn. (2) implies 
that: 
dU ЛІД Е и 


ту (0,1) 
dt Ye HNÉo и to (6) 


where 1,(0,t) is the stress of the elastic deformation at 

x = 0. For elastic energy accumulation in the deform- 

able domain, the inequality тє < 0 in Eqn. (6) occurs. 
If from Eqn. (6): 


_ 120.) _ 0,0 
УБ = . 


Then, Eqn. (6) can be rewritten as: 


dU Аб; U 
Pa T= (7) 
dt шо to 


furthermore Eqn. (7) can be written as: 


dU фр U 
ат (8) 
dt ct) ty 


where с = ишта. Equation (8), for a certain depend- 
ence of af) and according to the initial conditions, de- 
scribes the variation of the releasing energy stream time, 
for example during an earthquake. In the stationary case, 
U(t) satisfies the condition: 


dU 
—=0 9 
ЕР (9) 
which implies that: 
1 
Оз, == дет (10) 


The physical meaning of the coefficient 6, can be 
seen from Eqn. (10): i.e. it is the ratio of the releasing 
energy stream to the total stress. 
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Let the energy stream rise and damping be considered. 
At first, let there be no loading action applied to the to 
be investigated domain of rocks. Then assume a constant 
loading action (i.e. Е = 0, therefore т = const) at instant 
т = 0. For Eqn. (8), the initial condition is given as: 


U(0) =0 (11) 


Then the solution for Eqn. (8) can be written as: 


14| on : ) (12) 
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After the loading action is taken off (ie. т = 0), the 
energy stream damping corresponds to the initial condition: 


0(0) = 0, (13) 


According to Eqn. (13), from Eqn. (8) one obtains: 


U= vox- +) (14) 


0 


A comparison between Eqn. (12) and (14) shows 
that it is easier to measure tọ from the relaxation curve 
of the energy stream. 

Having obtained the energy stream relaxation curve, 
one can make the dependence In U/U,(t). It follows 
from Eqn. (14) that: 


t 


U 
te (15) 


i.e. the above dependence has the form of а straight li- 
ne. From the slope of this straight line, one obtains the 
viscoelastic deformation time tọ: 


= 
ty = -afa In z) (16) 
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The probability of viscoelastic deformation 


According to MELKONYAN (1997), elastic and viscous 
deformations probabilities are assumed as: 


Ке =—; Ry =— (17) 


where tog = (Ту = const» fov = (Dy = const, аге the elas- 
tic and viscous deformations times, respectively. The total 
probability of the viscoelastic deformation is assumed to 
be (e.g. MELKONYAN, 1997): 


t t 
TEE Еи (18) 
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Taking into account the deformation probability, the 
coefficient др by its physical meaning can be present- 
ed as the ratio of the elastic deformation probability to 
the total deformation probability, i.e.: 


ôg =-= (19) 


Combining Egns. (17), (18) and (19), one obtains: 
to 


ôg =— (20) 
tog 


According to Eqn. (12), for a loading duration At >> tọ 
one obtains: 


U=U,, (21) 


Substitution of Eqn. (21) into Eqn. (10) yields: 


U 
ar 
Е Е (22) 
ог, see Eqn. (20): 
U` of 
ЕЕЕ (23) 
т tor 


Viscoelastic deformation efficiency 


Let the ratio tọ/tọg be denoted by W; and called the 
efficiency of viscoelastic deformation 


W, =o (24) 


Substituting expression (18) into Eqn. (24), one obtains: 


toy 
г 25 
tog + toy ©?) 


Equation (23) implies that at At >> tọ, i.e. for a long 
duration of the loading, the ratio U/t depends on both 
the elastic and the viscous deformation times. 

At At << tọ, from Eqn. (12) one obtains: 


At 


U(t) =U gy, =. (26) 
0 


Substitution of Eqn. (10) into Eqn. (26), using Eqn. 
(20), implies that: 


= (27) 


It follows from Eqn. (27) that the dependence f9,(7) 
can be obtained from the dependence И/топ т. 
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Activation energy of viscoelastic deformation 


It is known that in solids as well as in liquids, the 
viscous flow is a diffusion phenomena. In this case the 
stress state leads to certain deflections of the direction 
of the diffusion components and the absolute quantity. 

The process is named viscous flow in liquids and 
pseudoviscous flow in solids (e.g. VERHOOGEN ef al., 
1974). The velocity of this process depends on the tem- 
perature. The process equation contains the Boltzmann’s 
constant k and the activation energy Е, (e.g. VERHO- 
OGEN et al., 1974): 


Еле Бо 
n р ЕТ (28) 


where А is a constant. 

However, according to Eqn. (4), from Eqn. (28), the 
temperature dependance of the viscoelastic deformation 
energy stream can be obtained: 


U~ e- 2, (29) 


For the case At << tọ, according to Eqns. (23) and (29), 
U depends on the temperature only through the elastic de- 
formation time ѓу. Therefore, the activation energy values 
obtained at Аг << tọ will be stipulated by the temperature 
dependence of the elastic deformation time: 


U~ an|- — | (30) 


where ЛЕ р is the activation energy of elastic deformation. 

The viscoelastic deformation energy stream ratios, 
obtained for the case of At >> К and At << tọ for the 
same value of 7, according to Eqns. (23) and (27), must 
be equal to the ratio АЛХ. Then, it can be assumed that: 


U(At >> to) to (31) 
U(At <<) At 

where U(At >> tọ) and U(At << tọ) are the elastic de- 

formation energy streams at an applied force duration 

of At >> fy and At << tọ, respectively. 

Taking into account that in rocks the deformation is, 
in principal, viscous (e.g. УЕВНООСЕМ et al., 1974), i.e. 
probabilities ratio др = R,/Rp << 1, and therefore 
to/ton << 1, see Eqn. (24), then from the temperature 
dependence of U(At >> tọ)/U(At << tọ), the activation 
energy of viscous deformation, ДЕ,у, is obtained i.e.: 


U~ an|- Е: | (32) 


Analysis of the dependence (29) shows that for a 
viscoelastic deformation energy stream at At >> tọ, tak- 
ing into account Eqn. (23), И ~ к, whence follows: 


AE AE_, -AE 
U~ ext = | = ex Е 14 (33) 


Thus, for the case of At >> tọ, the values of the acti- 
vation energy ДЕ, obtained correspond to the difference 
of the termal activation energies of the elastic and vis- 
cous deformations energy carriers: 


AE, = AE gg — AE gy (34) 


Conclusions 


The above-mentioned model, which considers the si- 
multaneous appearance of both elastic and viscous beha- 
viours of rocks, at the same time allows the elastic and 
viscous deformation times, as well as their activation 
energies to be separated. From this model many other 
characteristics of the viscoelastic deformation of rocks 
can be obtained. For example, from Eqn. (18), it can be 
seen that the energy stream relaxation time к may be 
increased both with respect to an increase of fg, and (у. 
Eqn. (24) shows that the viscoelastic deformation effi- 
ciency decreases with increasing elastic deformation time 
and increases with increasing viscous deformation time. 
From the dependence of the viscoelastic deformation effi- 
ciency W; оп f, one can determine which deformation 
time (elastic or viscous) is varied. If W, increases with 
fọ, then to, = const., if W; decreases, then Ку = const. 

The value of the activatin energy of the viscoelastic 
deformation AE, in the stationary case is connected with 
neither the elastic nor the viscous deformation activa- 
tion energies but is the result of the superposition of 
these energies, see Eqn. (34). 

The above methods for the definition of the charac- 
teristic time к, and the model allow the characteristics 
of the viscoelastic deformation to be determined, which 
could enable the viscoelastic deformation process to be 
controlled. 
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Резиме 


Мобилизација енергије при 
високоеластичним деформацијама стена 


Модел који је анализиран у овом раду приказује 
симултану појаву (у интермолекуларној размери) 
еластичног и вискозног понашања стена током њи- 
ховог вискоеластичног деформисања. На тај начин 
омогућено је дефинисање времена релаксације ty 
вискоеластичне деформационе енергије. Када се ово 
време одреди, онда могу да се раздвоје активиране 
енергије еластичних и вискозних деформација од 
температурне зависности вискоеластичног дефор- 
мационог енергетског тока. Пробабилистички при- 
ступ, примењен на вискоеластичне деформације 
стена, омогућава да се уведе појам тоталне вероват- 
ноће, која се састоји од вероватноће еластичних и 
вискозних деформација. 


На тај начин предложени модел разматра симул- 
тану појаву еластичног и вискозног понашања сте- 
на и у исто време омогућава да се раздвоји еластич- 
но и вискозно време деформисања као и њихове 
активиране енергије. Помоћу овог модела могу да 
се одреде и многе друге карактеристике вискоелас- 
тичног деформисања стена. На пример: енергетски 
ток времена релаксације f јача са повећањем како 
еластичног тако и вискозног времена деформације; 
вискоеластична деформациона ефикасност се сма- 
њује када се еластично деформационо време пове- 
ћава, а расте када се вискозно деформационо време 
повећава. На основу зависности између вискоелас- 
тичне деформационе ефикасности и (у, може да се 
одреди које се деформационо време мења (еластич- 
но или вискозно). Ако се деформациона ефикасност 
повећава са fọ, онда је еластично деформационо 
време константно, а ако се деформациона ефикас- 
ност смањује са tp, онда је вискозно деформационо 
време константно. Величина вискоеластичне де- 
формационе активиране енергије, у стационарном 
случају, није повезана ни са еластичним а ни са 
вискозним деформационим активираним енергија- 
ма, али је последица суперпозиција ових енергија. 


